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I. INTRODUCTION 



The wormhole concept was invented by Wheeler (see, e.g., [T]) to provide a mechanism for having charge 
without charge, since in a such a spacetime without a center, the field lines seen in one part of the Universe could 
thread the throat and reappear in other part. The idea, was further extended by Morris and Thorne to allow, not 
only field lines, but also observers to travel through the throat f2|. By having invoked an arbitrarily technologically 
advanced civilization, this work ,2j initiated the engineering of wormhole construction, theoretically, and the study 
of wormholes as topological bridges joining two different spacetimes has since then attracted extraordinary attention 
in this modern context (see, e.g., iSj and references therein). The main feature of traversable wormhole physics is 
the fact that the matter threading the wormhole throat should possess exotic properties, one of them being the 
violation of the null energy condition |3] . To provide the existence of wormholes one should cither include some 
hypothetical forms of matter into the model, or introduce interactions of a new type. Various models, admitting 
the required violation of the null energy condition, have been considered in the literature, among them thin shells 
in a cosmological constant background O [6j, scalar fields [7], wormhole solutions in semiclassical gravity 0, 
solutions in Brans-Dicke theory [9], wormholes on the brane 10], wormholes supported by matter with an exotic 
equation of state, namely, phantom energy [TT], the generalized Chaplygin gas 12 , tachyon matter [T3], nonlinear 
electrodynamics [TJ, and other cases. 

Now, when one considers nonminimal coupling of gravity with vector-type fields, such as Maxwell, Yang- 
Mills or Proca fields, new interesting possibilities appear. Nonminimal phenomena, i.e., phenomena induced by 
the interaction between curvature, or gravity, and other fields can be characterized on the one hand by unusual 
effective energy conditions, and on the other hand, allow one to exclude exotic substrates. In [15] an exact solution 
of the nonminimal Einstein- Yang-Mills model was obtained which demonstrates that an SU{2) symmetric gauge 
field nonminimally coupled to curvature can support a traversable wormhole. This is a nonminimal Wu-Yang 
magnetic wormhole, and the throat joins two asymptotically flat regions. The corresponding spacetime has no 
center and this model could be an illustration of Wheeler's idea about "charge without charge" [1]. Following this 
idea we now intend to consider an electrically charged object in the context of a spacetime without a center. A 
few wormhole models are known in which the electric charge is spread on the spherical shell |16j . or the throat 
is filled with some nonminimal and ghostlike scalar field [17j. Here, we find an exact solution of a nonminimal 
traversable wormhole, which contains neither electric charge on spherical shells, either thin or thick, nor scalar 
fields, nevertheless, possesses a static spherically symmetric electric field which is charged from the point of view 
of a distant observer. Thus, our goal is twofold. We present explicitly a nonminimal realization of Wheeler's idea 
about charge without charge, and we show that, if the world is somehow nonminimal in the coupling of gravity to 
electromagnetism, then wormhole appearance, or perhaps construction by an absurdly advanced civilization [5], is 
possible. 
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For this purpose we address a nonminimal Einstein-Maxwell theory, which has been earlier elaborated in detail 
in both linear (see, e.g., [I51[Tn]) and nonlinear (see [2D]) versions. The model linear in the spacetime curvature 
and quadratic in the Maxwell tensor |20| contains three nonminimal coupling constants gi, q2, and 93. These 
quantities have the dimensionality of area and characterize the cross terms in the Lagrangian linking the Maxwell 
field i^ifc and terms linear in the Ricci scalar i?, Ricci tensor Rik, and Riemann tensor Rikmm respectively. These 
parameters are a priori free ones but can acquire specific values in certain effective field theories. The first example 
of a calculation of the three coupling parameters was based on one-loop corrections to quantum electrodynamics 
in curved spacetime, a direct and nonphenomenological approach considered by Drummond and Hathrell |21j . In 
another instance, Buchdahl \7T and then Miiller-Hoissen [53^ obtained a nonminimal Einstein-Maxwell model from 
dimensional reduction of the Gauss-Bonnet action. This model contains one coupling parameter. 

Based on these works, and specially on [5D], solutions, which described nonminimal electrically charged objects 
with and without centers, characterized by two sets of relations for nonminimal coupling parameters, namely, 
<?! + 92 + 93 = 0, 2qi + 92 = 0, and gi -l- 92 = 0, — 0, were obtained However, none of these solutions 

can be used to construct traversable wormholes. In this paper we formulate a new nonminimal model, for which 
the coupling parameters satisfy the relations, 891 + 92 = 0, (73 = 0. Exact solutions of this model are shown to 
admit the existence of nonminimal traversable electrically charged wormholes. Comparing these new results with 
the ones obtained in |15| we would like to emphasize the following features. First, here we deal with an electric 
field instead of a magnetic one; second, in [TS] we used the basic conditions I2qi + 4^2 + 93 = 0, 93 ^ 0; third, the 
spacetime is now nonsymmetric, i.e., the throat joins one asymptotically flat region to another asymptotically de 
Sitter region. In [25! the authors have suggested the name "black universes" to nonsymmetric wormhole spacetime 
configurations for which the first asymptotics is fiat and the second one is of cosmological type. 

The paper is organized as follows: In Section |llj we quote the details of the three- parameter nonminimal 
Einstein-Maxwell model and formulate the corresponding key equations for the electric and gravitational fields. 
Furthermore, we introduce a new one-parameter nonminimal model, given by the conditions 3gi + 92 = 0, (73 = 0, 
derive a key (cubic) equation for the electric field E and obtain the metric functions cr, N in terms of E. In Section 
III, we discuss three solutions of the equation for the electric field. We focus on the Coulombian-type solution: This 
is the solution to which the others should be compared. In Section IV, we obtain exact solutions with a center. In 
Section V, we thoroughly discuss wormhole solutions. Section VI contains the conclusions. 



II. NONMINIMALLY EXTENDED EINSTEIN-MAXWELL THEORY 



A. Master equations 

Details of three-parameter nonminimal Einstein-Maxwell theory can be found in the papers [20), [24] ; here we 
extract the main elements of this model only. The action functional is of the form 



S'nmem — d xy/^ 



(1) 



Here, g — det(f/jfc) is the determinant of the metric tensor g^fc, R is the Ricci scalar, k, is the gravitational constant. 
The Latin indices run from to 3. The Maxwell tensor Fik is expressed, as usual, in terms of a potential four-vector 
Ak 

E,k = V.Afe - VkA, , (2) 
where the symbol denotes the covariant derivative. The tensor t^*'=™" is defined as follows (see [IS]): 

^ ikmn 9l jj / ira kn ^in km\ i / jjirn kn jjin km , jjkn im T:tkm in\ i „ jjikmn ( r>\ 

X = ^R{9 9 -99 ) + ^[R 9 - R 9 +R 9 -R 9 ) + q^R , (.•jj 

where R^'^ and are the Ricci and Riemann tensors, respectively, and qi, q2, 93 are the phenomenological 

parameters describing the nonminimal coupling of electromagnetic and gravitational fields. The variation of the 
action functional with respect to potential Ai yields 

VkH''' = , H''' = F'^ + x'^"''^Fran , (4) 

where is the nonminimal excitation tensor [26| . In a similar manner, the variation of the action functional 
with respect to the metric yields 



R^k^lR9^k^nTlf\ (5) 



The effective stress-energy tensor Tj-f^'' can be divided into four parts: 



7t"^=Tif^+9rTf^+92Tf^^+937f^^ (6) 
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The first term T- 
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ik — ^yik-^ mn-^ ^ in-^ k ' 



(7) 



is a stress-energy tensor of the pure electromagnetic field. The definitions of other three tensors are related to the 
corresponding coupling constants qi, q2, qs'- 



1 



1 



(8) 



T, 
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ik 



,9ik 
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-^9ik^ ^mn^ls—^^ \^ i ^knls ~r r ^inlsj^^^m^nY^i ^k ^ k ^ i \ 



(9) 
(10) 



One may check directly that the tensor T^^^^ satisfies the equation 'S/^T^'^^^ — 0. 



B. Spherically symmetric model with electric charge 

Below we consider the nonminimally extended Einstein-Maxwell equations, given in Eqs. (4)-(10), for the case 
of a static spherically symmetric spacetime, given by the metric 

ds" = a^Ndf -^-y^ {dO^ + sin' ^V) , (11) 

where N, a, and Y are functions of the radial variable r only satisfying the asymptotic flatness conditions when 
r — T' -l-oo 

a{r)^l, N{r)^l, Y{r) ^ r . (12) 

In this paper, asymptotic flatness means that we do not consider the case Y'{r) = identically. Let us assume also 
that the electromagnetic field inherits the static and spherical symmetries. Then the potential four-vector of the 
electric field Ai and the Maxwell tensor have the form 

A, = Aoir)S^ , F.,k = A',ir) {S[St ~ SlS[) , (13) 

where a prime denotes the derivative with respect to r. To characterize the electric field, it is useful to introduce 
a new scalar quantity E{r) as 

E^{r) ^ -^F,kF'^ (14) 

Then the electric field squared is (A'q)'^ from which one obtains in turn Fro = —cr{r)E{r). 
The Maxwell equations Q give only one nontrivial equation, namely, 

[Y^E{l + 2x''\o)]' = 0, (15) 

which can be integrated immediately to give 

1-^(1 +V°.o)' ^''^ 

where Q is a constant to be associated with an electric charge of the object. Supposing the spacetime to be 
asymptotically flat, i.e., RiMm{T ~^ +00) = 0, one can see that (16) yields asymptotically the Coulomb law 
E{r) — > Q/r^. The nonminimal excitation tensor H"^^ given in Q with x*^™" given by ^ has only one nonvanishing 
component H'^° . Let us introduce an electric excitation scalar as follows 

D = J-^H,kH^k = ^-HroH-^ = E{1 + 2x'"°,o) • (17) 



According to ( 16 1 the electric excitation takes a simple form 



D=§-,, (18) 



i.e., D also tends to zero at r ^ -foo. 
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C. Key equations 

1. Preliminaries 

Analogously to the static spherically symmetric model in minimal electrodynamics there are two independent 
Einstein equations only, e.g., 



- g; - At ( T,"^"''^ - T^^-^'A , GO = «;roO(°«) 



The first one gives the equation linking the function (7{r) with E{r) and Y{r) 
K{qi + q2 



a 



qs) [Y'^iEE'Y + YY"E^] + 2n{q2 + 2q:i)YY'EE' + Kggy'^i^^ _ yY" 
Y [n{qi +q2+ qs) [YEE' + Y'E^) - Y'] 



It is convenient to rewrite the second equation in ( 19 1 using ( 16 ) as 



Q 

Y^E 



1 



(19) 



(20) 



(21) 



For generic coupling constants this equation yields 

N'Y + 92 + qz){YEE' + Y'E^) ~Y']+N {2K(gi +q2+ q^) [Y^{EE'y + YY" E^] + 
+ 2KEE'YY'i2iq, +q2 + qs) + q2 + 2q^) + «;((gi + q2 + q^) + q3)E^Y'^ - Y'^ - 2YY"] + 

kE'^Y'^ 
2 



1 



- Kq^E'' - kQE = 0. 



(22) 



The three Eqs. (16), (20 1 and (22) form the key system of equations determining the four unknown functions cr. 



A^, E and Y . One of these functions is arbitrary. Generally, one has cr'(r) and N'{r) given in Eqs. (20) and (22) 
and put both into (16 1. This procedure yields 



2NE [K^iqi +q2+ qsfiSqi + q2)Y E^ {EE"YY' - EE'YY" - YY' E'^ + EE'Y'^) - 
- Kiqi + 92 + '?3)(3gi + 92 - q3)Y\EE"Y' - EE'Y" + Y'E'^) - 
- k('Zi + 92 + 93)(3<Zi + <Z2 - 'iq3)YY'^EE' + nq^iqi + q2 + q3)Y'^E^ - q^Y'^] - 
= {K^E^E'iq, +q2+ q3)[{2qi + q2)Y^ - 2q,Y{2qi + q2 + {qi + 92 + ^s))] 
~t^YQE^E'{q^ + 92 + 93)(4<Zi + 2q2 + {qi + 92 + 9,3)) - 
KE^E'iqi +q2+ q3){Y^ - 2Y{3qi + 92)) + nYQEE'{qi + 92 + 93) - K^Y'E''{qi + (72 + q3)qi{Y^ + 2(<Z2 + 93)) + 
+ K^QY'E^qi +q2 + g3)(<Zi - 92 - 93) - «1"'^'(92 + 293)(1^' - 291) + 2KQY'E^iq2 + 293) + 

+ Y'E{Y^ -2q3)-QY'} , (23) 



thus giving iV as a function of E, E' , E" , Y, Y' and F". In principle, putting N from Eq. (23 1 in Eq. (22) one 



can obtain a key equation linking E(r) and Y{r). This equation happens to be nonlinear and contains third order 
derivatives E'" and Y'" . Although there are no explicit exact solutions for arbitrary gi, 92, 93, one can estimate 
the behavior of the integral curves using a decomposition with respect to Y^^. For instance, when the model is 
asymptotically flat, the decompositions start as 



N 



E = 
a 

1 

yv2 



Q_ 

--Y' 



1 



14 

2M 
~Y~ 



493M _ K(g2 + 393)Q^ 
y3 

k(109i 



692 



5q3)Q^ 



2y2 



41-4 
k(49i 



392 + 293)Q' 



y4 



(24) 
(25) 
(26) 



where the constant M is an asymptotic mass of the object. 

There are two special cases arising when the left-hand side of Eq. ( 23 ) vanishes so that the function A'' cannot 
be found from this equation directly. The first one is the model with 91 + 92 + 93 = and 93 = 0, which has been 
studied in [53], the second relates to the model with 391 + 92 = and 93 = 0. Since the model with 91 + 92 + 93 = 
and 93 = has been studied in |24| we briefly mention that this one-parameter model can be characterized by the 
relationships 91 = —9, 92 
also vanishes, i.e.. 



9, 93 = 0, for some 9, and ( 23 ) is satisfled, when the right-hand side of this equation 
KqE^{Y^ + 2q) - 2KqQE^ - EY^ + Q = . (27) 
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The function a can be readily found as 



cr = y Cxp(-K(?£'2) 



(28) 



When the function Y coincides with the radius r, i.e., Y{r) = r, we obtain indeed the model which has been 
investigated in the paper 1 24j . On the other hand the model with 3gi + q2 — and = is also very interesting 
and will be studied in detail here. 



2. One-parameter model with 3qi + q2 ~ and qs = 

The one-parameter model with 3qi +(72 = and (73 = is characterized by the relationships qi — —q, q2 = 3q, 
(73 — 0, for some q. The susceptibility tensor is now proportional to the difference of the Riemann (Rikmn) and 
Weyl (Cikmn) tensors 



/Cikmn {^Rikrnn Cik'mri) 



Equation ( 23 ) is satisfied when 

- \Y' - 2Kq{YEE' + Y'E^)] [nqE^iY'^ + 6q) - AnqQE'^ - EY"^ + Q] = . 



(29) 



(30) 



Again, one has two variants to obtain exact solutions. The first one can be realized, when the first bracket vanishes, 
and we obtain 



Y' - 2Kq{YEE' + r'£;2) = ^ Y'^ = 



const 
1 - 2KqE^ 



(31) 



The minimal limit g — > relates to the model Y — > const instead of F — > r, and we do not consider here such a 
model. When the second bracket in (30) vanishes, we obtain the following cubic equation for the electric field 



KqE^{Y^ + 6q) - AnqQE^ - EY^ + = 0. 



(32) 



This equation differs from (27) by the numerical multipliers in front of q. Surprisingly, the equation for the function 

(33) 



a{r) [see Eq. (20)] can be now explicitly resolved in terms of E and Y as 

a = Y' - 2Kq (E^Y' + YEE') . 



In order to find the function N{r) we consider Eq. (22), which gives 

1 



N = 



g'^Y 



dra\\^ ' --"^^ 



aqE — nQE + const 



(34) 



Thus, any solution E{Y) of the cubic equation [32| for the electric field gives us a new exact solution of this 
nonminimal model, since (33) gives immediately the function a{r) and (34) gives the function N{r) in quadratures. 



3. Dimensionless quantities and equations 



introduce the following dimensionless quantities (see |Hj) 



In order to analyze Eqs. (|32|)-(|34|), and in particular the cubic equation (|32|) for the electric field, let us 

E 



Y 

rq rq Eq 



(35) 



In these terms the Coulombian branch of Eq . (|32[ ) corresponds to the solution with the asymptotic behavior Z ^ p ^ 
at p ^ +(X). Then the key equations (32)-([34|take the form 



a Z-\y^ + 3a) - 4a - Zy^ + 1 = , 



(36) 



dp 



2a ( Z'^'^ +yZ 



dp 



dZ 
dp 



(37) 



N = 



dpa (1 + Z'^y'^ + a Z'^ - 2Z) + const 



(38) 



The dimensionless parameter a is a guiding parameter for the key equation (36). When a = 0, then q = and the 
model is minimal. 
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4- Solution with a < 



Before we enter into the really interesting solutions, we mention the solutions with a < 0. The expression 
+ 3a, the coefficient of in the Eq. (36), can vanish, when a < 0; thus, the corresponding solution Z{y) is 
singular. The equation y = \/~3a describes a vertical asymptote. There exists one critical value of the guiding 
parameter a, namely, a — —16/9, for which the point characterizing by the condition — oo, belongs to this 
vertical asymptote. The typical behavior of electric field at a < is presented on Fig. [T] Anyway, we will consider 
regular solutions only, thus assuming that a > from now on. 




a) a < -16/9 




10- 



3 4 



5 6 



= -16/9 



c) -16/9 < a < 



FIG. 1: Plots of the rescaled electric field Z{y) = E/Eq as function of y = Y/fQ for a < 0. All curves have a vertical 
asymptote at y = y/—3a. (a) When — oo < a < —16/9, the curve crosses the asymptote twice, (b) When a = —16/9, 
the curve touches the asymptote at one point, (c) When —16/9 < a < 0, the curve does not cross the asymptote. For 
all cases the function Z{y) is infinite, the singularity point being situated at y = 0, when a = only. When a — >■ —0, the 
vertical asymptote coincides with the coordinate line y = and the right branch of the curve (c) converts into the standard 
Coulombian curve Z{y) = 1/y^. 



III. ELECTRIC FIELD: THE COULOMBIAN AND NON-COULOMBIAN SOLUTIONS 



Now we study the a > cases. The Coulombian solution for the electric field is the basic solution. Indeed, 
there are three branches of solutions to Eq. (36 1 with a > 0. One of them Z^{y,a) can be identified as the 
Coulombian branch, the other two can be indicated as the first and second non-Coulombian solutions, Z_(y,a) 



and Z+{y,a), respectively. Let us study first the Coulombian solution. From Eq. (36) this solution is 




12a 



243a 
^56" 



( 2y« 



^ - 1 
9a ^ 



3/2 



> . 



(39) 



Clearly, Z^,{y,a) — ?> y~'^, when y — ?> oo for arbitrary a, thus, this solution can indeed be referred to as the 
Coulombian branch. This solution remains real when the modulus of the arccosine argument does not exceed 
unity, i.e., y belongs to the interval given by the inequality 



243a 



2f_ 
81a2 



r 

9a 



- 1 



< 1 



y 

12a 



3/2 



(40) 



In particular, when a < 256/243, then this interval is < 2/ < oo- At y = this solution is not singular, since 



^*(0,a) 



243a 
^56" 



< a < 



256 
243 



(41) 



The function Z^,{0, a) as a function of the guiding parameter a has a minimum at a = 8/9, the minimal value being 
2^(0)inin = 3/4. The behavior of the Coulombian branch of the solution for E{r) can be illustrated in a figure. In 
Fig. 2^-e we plot the rescaled electric field Z{y) as a function of the rescaled radius y, for different values of a. The 
Coulombian branch Z^{y,a) is represented by the middle line. 

The second solution Z_(y,a) of the cubic equation (36) can be represented in the form 
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d) a = 256/243 e) a > 256/243 



FIG. 2: Plots of the rescaled electric field Z{y) displaying the behavior of the Coulombian branch and the two non- 
Coulombian branches for a > 0. The interval < a < +oo is split into three regions with respect to two specific values of 
a, a = 1 and a = 256/243. The curve Z{y) is also drawn for the cases a = 1 and a — 256/243. Dashed lines indicate upper 
(Z = 1/^/a) and lower {Z = —1/^/a) horizontal asymptotes, (a) When < a < 1, Z+{y,a) > (b) When a = 1, the 

second non-Coulombian branch degenerates into a straight line Z+{y, 1) = I. (c) When 1 < a < 256/243, Z+(y,a) < 1/^/a. 
(d) When a = 256/243, the second non-Coulombian branch has a common point with the Coulombian one at y = 0. (e) 
When a > 256/243, these branches have a common point at j/ > 0. 



At ?/ — >■ oo this curve tends to the horizontal asymptote Z = — 1/y^, and in this sense it is non-Coulonibian, we call 
it the first non-Coulombian solution. The solution is real for arbitrary (positive) value of the guiding parameter a, 
see the lower curves of the Fig. [2k-e. At the point y ~ Q, the center, this solution is also nonsingular 



Z_(0) 
Z_(0) 



1 

'4 V a 
-27/64 
1 [3 



cosh 



1 



arccosh 



/243a 
^56" 



/243a 
^56" 



a < 256/243, 
a = 256/243, 
a > 256/243. 



(43) 



The third solution of the cubic equation ( 36 ) is of the form 




12a 



- arcsm 



243a 



1 - 



y_ 

9a 



81a2 



1 



12a 



3/2 



(44) 



For y ^ oo this curve tends to the horizontal asymptote Z — l/ya, thus the solution is also non-Coulombian, it 
is the second non-Coulombian solution. Nevertheless, now the curve behaves analogously to the Coulombian one, 
for instance, the solution is real for < y < +oo when a < 256/243, see the upper curves of the Figs. ^-e. The 
value at the center is 



Z+{0,a) = 



1 



1 



■ arcsm 



/243a 
^56" 



a < 



256 
243 



(45) 



Clearly, Z*(0) — Z+{0) — 27/32, when a — |||, i.e., the curves and Z+ contact at ?/ = at this value of the 
guiding parameter a, see Fig.[2]3. When a > 256/243 the curves Z^ and Z^ contact at y > 0, the inverse function 
y{Z) being smooth in the vicinity of this point. Thus, when a > ||| one has that ymin > 0, and there is the 
possibility of finding a throat for a traversable wormhole. 
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Let us also mention that Z_|_(0,a) > Z_|_((X),a) when a 



< 1, and Z_|_(0,a) < Z+(oo,a) when a > 1 (see 



Figs. [2]a, c). In this connection we would like to make the following remarks: The cubic equation (36) differs from 
the analogous equation discussed in [21] [see (27)] by the coefficients in front of and . This leads to the 
following differences: ('ij In contrast to [21] the curves and Z+ do not intersect at a = 1. (ii) There are no 
bubbles in the vicinity of the center in the framework of this new model. (Hi) The branches Z^ and can contact 
at y = 0, it is possible, when a = 256/243. 



IV. EXACT SOLUTIONS: SOLUTIONS WITH A CENTER 



A. Behavior of the function cr(p) 



Based on the study of the cubic equation for the electric field, we consider the solutions for the metric functions 
(j{r) and N{r). We focus here on the solutions with center, i.e., on the cases, when the point y = is accessible. 
The solution with center for the nonminimal electric field with Coulombian asymptote exists, when < a < 256/243 
(see Figs.[2^-d). Of course, solutions with center with non- Coulombian asymptotes exist for other a. 

Let us consider the behavior of the function a{p) for the case < a < 256/243. We can put now Y — r and, 
respectively, y = p ^ , keeping in mind the asymptotic behavior (12). Then the formula (37) takes the form 



a^l-2a[Z^ + pZ 



According to Fig. ji^-d, Z{0) and ^^^^ are finite and the value of the function al^p) at the origin reads 

cr(0,a) = l-2aZ^{0,a). 



Substituting (41 ) into (47) we obtain for < a < ||| 



cr(0,a) = 1 




(46) 



(47) 



(48) 



0, 



The value <7(0, a) as a function of the guiding parameter a is steadily decreasing with increasing a. When c 
this function tends to 1/2, when a — 256/243 then ct(0, a) = —1/2, and when a = 8/9 the function (48) takes the 



value zero. When 8/9 < a < 256/243, the function (7(0, a) is nonpositive. Thus taking into account the asymptotic 
condition a{oo,a) = 1 > there is at least one root Ps > the equation a{p) = 0. When a < 8/9, the value 
(t(0, a) is positive. 



B. Exact solution for a — 1 



When a — 1 Eq. (36) admits the solution Z — 1, which it is one of the curves from the family Z^{y,a), 



namely, Z^{y, 1). Thus, in this case the cubic equation (36 1 splits, so that the Coulombian solution and the 
first non-Coulombian solution Z_ can be obtained as sol 



utions of the corresponding quadratic equations. This 



procedure yields. 



Z*iy,l) 



\/l3 + V 



r 



2(3 + y2) 



Z.(0,1) 



13 



(49) 



Z-{y,i) = 



Vl3 + 2t/+y4^1_,^2 

2(3 + y2) 



Z-{0,1) = 



1 - V13 
6 



(50) 



We now study the behavior of the functions cr{p) and N{p) for this model with a — 1 which is of interest for 
illustrative purposes. In this case the Coulombian branch of the electric field is described by (49), thus the metric 
function a{p) is of the form 



(6 + 34p2)^13 + 2p2 



p4 - 39 + 82p^ + lOp" + lOpt' 
(3 + p2)3^i3 + 2p2T^ 



(51) 



The plot of the function a{p) is represented in Fig. [s] At the origin cr(0, 1) = — (\/l3 — 2)/9 ~ —0.178 is negative. 
This function has only one root ps = 0.289. At p = 1.706 it reaches a maximum, (Tmax = 1.067, and then tends to 
one asymptotically. 
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d) m = mo = 0.7924 



e) m > mo 



FIG. 4: Plots of the function N{p) for a = 1 for various values of the rescaled asymptotic mass m = M/rq. We split the 
interval < m < +oo into three regions with respect to two specific values of mass: m = m* = 0.7913 and m = mo = 0.7924. 
All curves except (b) possess a vertical asymptote a,t p — ps = 0.289. When m = mo [plot (d)] A''(0) is finite, contrary to 
the other cases. 



The constant in (38) is connected with the asymptotic mass M defined by M=^ lini [r(l— goo)]; where goo 



a N. Taking into account this condition one obtains 



500 = 1 - 



2m - J{p) 



N = 



1 - 



2m " J(p) 



J(p) ^ J dp 

p 



2(3 + p- 



""-^ (l9 + 23p2 + + p6 - (5 + 2p2 + p4)^i3 + 2p2 + p4^ 



(52) 



where m = M/rq. The value of the integral J(p) at p = is finite and equal to J(0) = 1.5848. Depending on the 
value of mass m there are three variants of the behavior of the functions gooip) and N{p) in the vicinity of p_= 0. 
They are: (i) m < mo = ^^(0) = 0.7924. The function goo tends to positive infinity, when p — >• [see Fig. ^a-c]. 
(ii) m = mo- The function goo takes the negative finite value (?oo(0) = —0.0096 [see Fig. p^d]. (in) m > mo- The 
function goo tends to negative infinity, when p — >■ [see Fig. [5]e] . Plots of the function N\p) and of the function 

500(P) " . — rn .rn 



for a = 1 are given in Figs. [4] and [sj respectively. 



In addition to the special value of the mass m = mg, there exists another special value m* — 0.7913, which 
can be obtained as follows. The function a{p) takes zero value p — ps (see Fig.p|, and N{p) — goo/o"^ generally 
becomes infinite at this point. Nevertheless, the possibility exists to make N{ps) finite, when ps is a double root of 
the function goo- This situation takes place, when m = m*, see Fig. |4]b for N{p) and Fig.[5]b for goo{p)- However, 
for both cases the Kretschmann scalar K. = RikmnR^'^™'"' is infinite at this point, K. oc {p — ps)^^ when m m* and 
JC (X {p — Ps)~^ when m = m* . Hence, the sphere with the radius p = ps is a singularity. In addition there exists 
a second singularity at the origin. It is of a Schwarzschild type, when m ^ mo, and becomes a conic singularity. 
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d) m = mo = 0.7924 




e) m > mo 



FIG. 5: Plots of the function goo{p) ~ o-^N for a = 1 for various values of the rescaled asymptotic mass m — M/rq. We 
spht the interval < m < +oo into three regions with respect to two specific values of mass: m = m* = 0.7913 and 
m = mo = 0.7924. (a) When m < m*, the function goo is positive and tends to infinity at origin, (b) When m — m* , the 
curve touches the horizontal axis at p — ps = 0.289. (c) When m* < m < mo, the curve intersects the horizontal axis at two 
different points; the function goo is negative between these points, (d) When m = mo, the curve intersects the horizontal 
axis at one point; the function goo tends to a finite negative value at origin, (e) When m > mo, the function goo tends to 
negative infinity at origin. 



when m = mo, since the functions a and N are finite at p = 0, but N{Q) ^ 1. When m > m* these singularities 
are hidden inside an event horizon, while for the case m < m* they are naked, see Figs. [4] and [Sj This behavior 
difi^ers from the situation described in ]T7\ for the magnetic monopole with the same set of nonminimal coupling 
constants. 



V. EXACT SOLUTIONS: WORMHOLES (SOLUTIONS WITHOUT A CENTER) 

A. The expression for Y 

According to Fig. [2^, when a > 256/243, the value F = is inaccessible for the Coulombian branch of the 
electric field. It means that Y{r) reaches a minimal value y,„in(a), for which the Coulombian branch smoothly 
joins with the second non- Coulombian one. Therefore we are in the presence of a wormhole. The corresponding 
surface of minimal area AttY^^^^ is its throat. 

In order to describe this object let us introduce a new radial coordinate r, which covers the entire range 
(— oo; +oo). Without loss of generality, we can take as usual the throat to occur at r = 

y'(o) = o, y"(o)>o. (53) 

Also we assume that positive values of r are associated with the Coulombian branch of the electric field, while the 
range (—00; 0) corresponds to the non-Coulombian one we are studying. We have to supplement the asymptotic 



form (12) for the function Y{r) with a condition at negative infinity 

lim Y{r) = +00 , (54) 



which correlates with Fig. [2^. 
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B. Radius of the throat 



Further analytical progress is possible if we fix an explicit expression for this function, satisfying both the 
asymptotic condition ( 54 1 and the throat conditions ( 53 1 . The simplest function of such a type is 



Y{r) = + r|62 , 



or, after rescaling ( 35 ) 



v{p) 



Y 



(55) 



(56) 



where h is the dimensionless radius of the throat. The value of h can be found as follows: In Eq. (39) the argument 
of the function arc cosine takes the maximum value +1 corresponding to y cx). The minimal value —1 relates to 
V = J/min = Thus the value of the throat's radius satisfies the equation 



61 



3a6^ + ^a&^ + ya^&^ - y (243a - 256) = . (57) 

Let us mention that the argument of the arc sine in Eq. (44) takes the value +1, and Z*(6,a) = Z^[b,a), as it 
should be at the junction point. Equation (57) is quartic with respect to b"^. For a > 256/243 this equation admits 



only one positive solution. The value of the throat's radius as a function of the nonminimal parameter a is presented 
in Fig. ^ When a — 256/243, b vanishes. When a » 1, this function behaves as 5(a) — 3 (|)^^^ 



5_ 
12 



0(1). 




FIG. 6: Plot of the dimensionless radius of the wormhole throat 6 as a function of the parameter a. It starts at a = 256/243, 
i.e., if the parameter a is less than this value a throat does not exist. 



For illustration, it is convenient to consider the value a = 4, since 6(4) = 2 and Z^^(6, 4) = 4 [see (p39|]. 
Substituting y{p) taken from ( [56| into the cubic equation (36) we obtain a plot of Z(p, 4) presented in Fig. [t] This 
curve is nonsymmctric: at p — >■ +cxd the function Z{p) tends to zero (Coulombian branch), at p — > — oo it tends to 
1/2 (second non-Coulombian branch). 



C. Asymptotic behavior 



When p — )■ +00 the solution for the electric field behaves according to the Coulombian law. In this limit the 
functions a{p) and N{p) go to one [according to ( 12 ) , ( 25 ) , and ( 26 )] , and we deal with the traditional asymptotically 
flat Minkowski spacetime. When p — oo the electric field Z(p) tends to the asymptotic value l/\/a 7^ 0. In this 



case a{p) tends to one as well, since in (37) ^ — J' —1, while N{p) becomes infinite 



Nip) 



3a 



(58) 



Thus, for p — > —00 the metric of this spacetime is asymptotically de Sitter, and the Riemann tensor in this limit 
may be decomposed as 

,.2 



R 



= -K {5l,5i - S\,5'^:) + o(r-i) , K = 



'Q 



1 

3a 6q 

The effective A-term induced by the electric field nonminimally coupled to gravity is thus equal to Acff 
l/{2q). In contrast to the electric strength E, the electric excitation D [see Eq. (17)], being equal to 

2 



(59) 
= 

KQip^ + 62) ' (^^^ 

tends to zero both for p — +00 and p — > —00. Similar situations, where the excitation tensor was equal to zero, 
while the electric, magnetic, or other fields did not vanish, were considered earlier in |28) . 



D = 



Y2 
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FIG. 7: Plot of the rescaled electric field Z{p) for a — 4. The region p > is the Coulombian branch [Z{+<x) — 0], the 
region of negative values of p is the second non-Coulombian one [Z{—oo) — 1/2]. The junction point p = is the point of 
inflection on this curve. 



D. Horizons and singularities 



If we intend to consider traversable wormholes only, we have to require the spacetime metric to possess no 
singularities. This means that the function cr(p) should not vanish anywhere, because, when cr(p) — 0, we have 
a singularity, as mentioned previously. Therefore, for the model considered in [24] with qi = —q, 92 = 9, 93 = 
traversable wormholes are absent, since in the throat ct = (see Eq. (28)). In contrast to [23], for our case (gi = — g, 
92 = 3g, 93 = 0) we have, 

a(0) = ~2abZ[0) > , (61) 

dp 

because ^ < at p = (see, e.g.. Fig. ItI. However a singularity can be located at p 7^ 0. In order to find the 
range for the parameter a for which the function a{p) is positive everywhere, let us consider the situation where 
the function a{p) given by (37) with Z{p) from (36) has a double root [in this case the curve of a{p) touches the 
horizontal axis]. Numerical calculations show that this is possible for two values of the parameter a: a = ai « 1.072 
and a = 02 ~ 39.380. When 256/243 < a < ai, the equation cr(p, a) ~ has two positive roots; when a > a2 
the two corresponding roots are negative. When ai < a < 02, the equation a{p,a) — has no real roots, i.e., 
cr > for any value of r. Figure |8] illustrates these features for a typical value a = 4, which is in the interval 
1.072 < a < 39.380. 




FIG. 8: Plot of the metric function (j{p) for a = 4. This function is positive and tends to unity at both infinities. The curve 
is nonsymmetric with respect to the vertical axis. 



Now, since the wormhole spacetime metric is an asymptotically de Sitter one, it possesses at least one horizon 
(the discussion of some aspects of wormhole physics in a cosmological context can be found, e.g., in [5] [6]). Thus, 
our wormhole configuration may be considered as traversable, if, in addition to the requirement of the absence of 
singularities, we suppose that horizons of the metric are out of a traveler way through the throat. For the metric 
function N{p) it means that in the vicinity of the throat N{p) is positive. To illustrate the behavior of the functions 
gm{p) and N{p) let us consider again the model with a — A, see Fig. [9j When the dimensionless mass parameter 
m is less than a critical value m (being equal to 1.0081 for a = 4), an event horizon is located behind the throat, 
hence in this case the wormhole throat is traversable. When m > to, a horizon is situated in front of the throat, 
and a traveler cannot go through, it is a nontraversable wormhole. 



E. Comparisons 



It is of interest to briefly compare our electric wormhole solutions with the magnetic nonminimal Wu-Yang 
wormhole solutions studied in |15j . The magnetic nonminimal traversable wormhole joins two asymptotically 
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FIG. 9: Plots of (upper row) and N{p) (lower row) for the value a = 4, wheu there is only one horizon. We split the 
interval < m < +oo into two regions with respect to a specific value of the rescaled mass m = 1.0081. (a) When m < m, 
the curves cross the horizontal axis at negative values of p, so that the horizon is located behind the throat, (b) When 
m = rh, then 300(0) = and iV(0) = 0, thus, the horizon is located on the throat, (c) When m > fh, the horizon is located 
in front of the throat. Plot (a) illustrates the case in which the horizon belongs to an asymptotically de Sitter region of the 
combined spacetime, it is a standard cosmological apparent horizon. In this sense such a wormhole is traversable. 



Minkowskian flat regions and the magnetic gauge field is symmetric with respect to wormhole throat, i.e., this field 
asymptotically vanishes in each region of the combined spacetime. On the other hand, the electric nonminimal 
wormhole joins an asymptotically flat region and an asymptotically de Sitter region, for which the corresponding 
electric field tends to a constant value. This introduces a new efi^ective, nonminimally induced, cosmological 
constant A^g = l/{2q). 



VI. CONCLUSIONS 



We have presented new exact solutions of a nonminimal Einstein-Maxwell model with coupling parameters 
satisfying the relations qi = —q, 92 = 3(7, qz = 0, with q arbitrary, which describe spherically symmetric electrically 
charged objects with and without center. The electric field of the objects is characterized by the following features: 
(i) it satisfies a cubic key equation and splits into three branches, one of them has Coulombian asymptotics, 
E{r) Q/r^; (ii) when the guiding nonminimal dimensionless parameter a = Aq/nQ'^ is positive, the electric 
field is described by a smooth function finite everywhere; (iii) when a > 256/243 the Coulombian branch of the 
electric field conjugates with the non- Coulombian branch with constant asymptotics E — the position of 

the junction point depends on the parameter a and coincides with geometrical center for a — 256/243. 

There are wormhole solutions. Indeed, when the nonminimal guiding parameter a is in the interval 1.072 < a < 
39.380 there are no spacetime singularities and the nonminimal field configuration without center is a wormhole. 
This is an explicit example of a nonminimal traversable electrically charged wormhole joining two regions of 
spacetime, an asymptotically fiat region and an asymptotically de Sitter region. The spacetime of this wormhole 
has at least one horizon, which, depending on the value of the rescaled asymptotic mass m, can be situated in front 
of the throat, just on the throat and behind it. When m, docs not exceed some critical value (e.g., m < 1.0081 for 
a = 4), the horizon is located behind the wormhole throat, i.e., in the asymptotically de Sitter region of the combined 
spacetime. Such a configuration is thus a traversable wormhole supported by an electric field nonminimally coupled 
to gravity. In this manner we have presented explicitly a nonminimal realization of Wheeler's id(^a about charge 
without charge, and showed that, if the world is nonminimal in the coupling of gravity to electromagnetism, then 
wormhole appearance , or perhaps construction by an absurdly advanced civilization, is possible. 

Perhaps, a next natural step is to consider nonminimal Einstein-Maxwell models with an a priori cosmological 
constant A. It is expect that one can also find solutions of our nonminimal Einstein-Maxwell model. These solution 
would describe electric or magnetic wormholes joining different combinations of asymptotically de Sitter, anti de 
Sitter, and Minkowski regions. For instance, an anti de Sitter-Minkowski electric nonminimal wormhole is expected 
to be free of horizons. 
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